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Abstract 

We have studied the leptonic decay Df — ► t + v t , via the decay channel r + — > e + f e P T , using a 
sample of tagged -D+ decays collected near the D*^Df peak production energy in e + e~ collisions 
with the CLEO-c detector. We obtain B{D+ -> r+i/ r ) = (5.30 ± 0.47 ± 0.22)% and determine the 
decay constant fo s = (252.5 ± 11.1 ±5.2) MeV, where the first uncertainties are statistical and the 
second are systematic. 



* Deceased 
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I. INTRODUCTION 



The leptonic decays of a charged pseudoscalar meson P + are processes of the type P + — > 
£ + z/£, where £ — e, fi, or r. Because no strong interactions are present in the leptonic 
final state £ + vi, such decays provide a clean way to probe the complex, strong interactions 
that bind the quark and antiquark within the initial-state meson. In these decays, strong 
interaction effects can be parametrized by a single quantity, fp, the pseudoscalar meson 
decay constant. The leptonic decay rate can be measured by experiment, and the decay 
constant can be determined by the equation (ignoring radiative corrections) 

T(P Q , tv,) = G ^f f K P m} (l " ^|) 2 , (1) 

where Gp is the Fermi coupling constant, Vq q is the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix p], 0] element, mp is the mass of the meson, and is the mass of the charged 
lepton. The quantity fp describes the amplitude for the Q and g-quarks within the P to 
have zero separation, a condition necessary for them to annihilate into the virtual W + boson 
that produces the £ + vi pair. 

The experimental determination of decay constants is one of the most important tests of 
calculations involving nonperturbative QCD. Such calculations have been performed using 
various models [3J| or using lattice QCD (LQCD). The latter is now generally considered to 
be the most reliable way to calculate the quantity. 

Knowledge of decay constants is important for describing several key processes, such as 
B — B mixing, which depends on .fp, a quantity that is also predicted by LQCD calculations. 
Experimental determination 0, [1| of fp with the leptonic decay of a B + meson is, however, 
very limited as the rate is highly suppressed due to the smallness of the magnitude of the 
relevant CKM matrix element V u b- The charm mesons, D + and D+, are better instruments to 
study the leptonic decays of heavy mesons since these decays are either less CKM suppressed 
or favored, i.e., T(D + -> £+v t ) oc \V cd \ 2 « (0.23) 2 and V{D+ -> £ + v t ) oc \V CS \ 2 « (0.97) 2 are 
much larger than T(B + — > t^vi) oc jV^I 2 ~ (0.004) 2 . Thus, the decay constants fp and fp g 
determined from charm meson decays can be used to test and validate the necessary LQCD 
calculations applicable to the B- meson sector. 

Among the leptonic decays in the charm-quark sector, — > £ + vi decays are more 
accessible since they are CKM favored. Furthermore, the large mass of the r lepton removes 
the helicity suppression that is present in the decays to lighter leptons. The existence of 
multiple neutrinos in the final state, however, makes measurement of this decay challenging. 

Physics beyond the standard model (SM) might also affect leptonic decays of charmed 
mesons. Depending on the non-SM features, the ratio of T(D + — > £ + ui)/T(D^ — > £ + Vi) 
could be affected [a], as could the ratio 0, @| — > t + v t ) /T(D+ — ► /i + z/ M ). Any of the 
individual widths might be increased or decreased. There is an indication of a discrepancy 
between the experimental determinations [3[ of fp s and the most recent precision LQCD cal- 
culation This disagreement is particularly puzzling since the CLEO-c determination [l(| 
of fp agrees well with the LQCD calculation [9| of that quantity. Some conjecture that 
this discrepancy may be explained by a charged Higgs boson or a leptoquark. 

In this article, we report an improved measurement of the absolute branching fraction of 
the leptonic decay — > t + u t (charge-conjugate modes are implied), with r + — > e + z/ e z/ T , 
from which we determine the decay constant fp g . 
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II. DATA AND THE CLEO-c DETECTOR 



We use a data sample of e + e~ — > D^DJ events provided by the Cornell Electron Storage 
Ring (CESR) and collected by the CLEO-c detector at the center-of-mass (CM) energy 
4170 MeV, near D^DJ peak production [l^]. The data sample consists of an integrated 



luminosity of 602 pb _1 containing 5.5 x 10 5 D^DJ pairs. We have previously reported [1^ 
14j measurements of — > /U + z/ M and — > t + v t with a subsample of these data. A 
companion article [TEJ reports measurements of fjj s from Df — > /U + z/ M and Df — > r + i/ T , with 



7T 



u T , using essentially the same data sample as the one used in this measurement. 



The CLEO-c detector [16l . 1 171 1 181 . 1 1 91 ] is a general-purpose solenoidal detector with four 



concentric components utilized in this measurement: a small-radius six-layer stereo wire drift 
chamber, a 47-layer main drift chamber, a Ring-Imaging Cherenkov (RICH) detector, and 
an electromagnetic calorimeter consisting of 7800 CsI(Tl) crystals. The two drift chambers 
operate in a 1.0 T magnetic field and provide charged particle tracking in a solid angle of 
93% of Att. The chambers achieve a momentum resolution of ~ 0.6% at p — 1 GeV/c. The 
main drift chamber also provides specific-ionization (dE/dx) measurements that discriminate 
between charged pions and kaons. The RICH detector covers approximately 80% of 4-7T 
and provides additional separation of pions and kaons at high momentum. The photon 
energy resolution of the calorimeter is 2.2% at E 7 — 1 GeV and 5% at 100 MeV. Electron 
identification is based on a likelihood variable that combines the information from the RICH 
detector, dE/dx, and the ratio of electromagnetic shower energy to track momentum (E/p). 

We use a GEANT-based 20j Monte Carlo (MC) simulation program to study efficiency 
of signal-event selection and background processes. Physics events are generated by evt- 

GEN 



21]. tuned with much improved knowledge of charm decays [22|, 23[, and final-state 



radiation (FSR) is modeled by the photos [24 
radiation (ISR) is based on cross sections for D 



program. The modeling of initial-state 
Df production at lower energies obtained 



from the CLEO-c energy scan [12j near the CM energy where we collect the sample. 



III. ANALYSIS METHOD 

The presence of two Df mesons in a D^DJ event allows us to define a single-tag (ST) 
sample in which a D J is reconstructed in a hadronic decay mode and a further double-tagged 
(DT) subsample in which an additional is required as a signature of decay, the e ± 
being the daughter of the r ± . The D~ reconstructed in the ST sample can be either primary 
or secondary from D* s ~ — > (or D*~ — > tt°D~). The ST yield can be expressed as 

n ST = 2NB ST e ST , (2) 

where N is the produced number of D^DJ pairs, £>st is the branching fraction of hadronic 
modes used in the ST sample, and esT is the ST efficiency. The ust counts the candidates, 
not events, and the factor of 2 comes from the sum of and D~ tags. 

Our double-tag (DT) sample is formed from events with only a single charged track, 
identified as an e + , in addition to a ST. The yield can be expressed as 

n D T = 2A^i3 S T^LeDT, (3) 

where £>l is the leptonic decay branching fraction, including the subbranching fraction of 
t + — > e + v e v T decay, £dt is the efficiency of finding the ST and the leptonic decay in the 
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same event. From the ST and DT yields we can obtain an absolute branching fraction of 
the leptonic decay £>l, without needing to know the integrated luminosity or the produced 
number of D^DJ pairs, 

n DT e ST n DT /e 

Oh = = , (4J 

"-st £dt ^st 

where e (= gdt / £st) is the effective signal efficiency Because of the large solid angle ac- 
ceptance with high segmentation of the CLEO-c detector and the low multiplicity of the 
events with which we are concerned, €dt ~ ^st^l, where cl is the leptonic decay efficiency. 
Hence, the ratio €bt/^st is insensitive to most systematic effects associated with the ST, 
and the signal branching fraction £>l obtained using this procedure is nearly independent of 
the efficiency of the tagging mode. 



A. Event and tag selection 

To minimize systematic uncertainties, we tag using three two-body hadronic decay modes 
with only charged particles in the final state. The three ST modes 1 are D~ — > 4>tt~ , D~ — > 
K~K*°, and D~ — > K S K~. Using these tag modes also helps to reduce the tag bias which 
would be caused by the correlation between the tag side and the signal side reconstruction 
if tag modes with high multiplicity and large background were used. The effect of the tag 
bias 6 t ag can be expressed in terms of the signal efficiency e defined by 

e DT e D Te'sT £ l4t e ST c h /r\ 
e = = = — = e iAag, (5) 

6 S T ^ST e ST &3T esT 

where e' ST is the ST efficiency when the recoiling system is the signal leptonic decay with 
single e in the other side of the tag. As the general ST efficiency esT; when the recoiling 
system is any possible D s decays, will be lower than the e' ST , sizable tag bias could be 
introduced if the multiplicity of the tag mode were high, or the tag mode were to include 
neutral particles in the final state. As shown in Sec. HVl this effect is negligible in our chosen 
clean tag modes. 

The Kg — > 7r + 7r~ decay is reconstructed by combining oppositely charged tracks that 
originate from a common vertex and that have an invariant mass within ±12 MeV of the 
nominal mass (||. We require the resonance decay to satisfy the following mass windows 
around the nominal masses Q: ^ K + K~ (±10 MeV) and K*° -> K + n- (±75 MeV). 
We require the momenta of charged particles to be 100 MeV or greater to suppress the slow 
pion background from D*D* decays (through D* — > ttD). We identify a ST by using the 
invariant mass of the tag M(D S ) and recoil mass against the tag M recoi i(D s ). The recoil 
mass is defined as 

M recon {D s ) = ^J{E ee -E Ds f- |p e e-pDj 2 , (6) 

where (E ee , p ee ) is the net four-momentum of the e + e~ beam, taking the finite beam crossing 
angle into account; (Ex) s ,pn ) is the four-momentum of the tag, with Ed s computed from 
Po s an d the nominal mass [3| of the D s meson. We require the recoil mass to be within 55 



1 The notations D~ — > <frir~ and D~ — * K~K*° are shorthand labels for D~ — > K~K + ir~ events within 
mass windows (described below) of the (f> peak in M(K + K~) and the K*° peak in M(K + tt~), respectively. 
No attempt is made to separate these resonance components in the K + K~n + Dalitz plot. 
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FIG. 1: The mass difference AM(D S ) = M(D S ) — m£> g distributions in each tag mode. We fit the 
AM(D S ) distribution (open circle) to the sum (solid curve) of signal function (double Gaussian) 
plus background function (second-degree Chebyshev polynomial, dashed curve). 



TABLE I: Summary of single-tag (ST) yields, where n| T is the yield in the ST mass signal region, 
ng T is the yield in the sideband region, s is the sideband scaling factor, and nsT is the scaled 
sideband-subtracted yield. 



Tag mode 


n ST 


n ST 


s 


re S T 


D- - 


10459 


807 


0.980 


9668.1 ± 106.1 


D- K~K*° 


18319 


7381 


1.000 


10938.0 ± 160.3 


D 7 K-K° s 


7135 


1409 


0.999 


5727.8 ± 92.4 


Total 








26333.9 ±213.3 



MeV of the D* mass [3J. This loose window allows both primary and secondary D s tags to 
be selected. 

To estimate the backgrounds in our ST and DT yields from the wrong tag combinations 
(incorrect combinations that, by chance, lie within the AM(D S ) signal region), we use the 
tag invariant mass sidebands. We define the signal region as —20 MeV < AM(D S ) < +20 
MeV, and the sideband regions as -55 MeV < AM(D S ) < -35 MeV or +35 MeV < 
AM(D S ) < +55 MeV, where AM(D S ) = M{D S ) - m Ds is the difference between the tag 
mass and the nominal mass. We fit the ST AM(D S ) distributions to the sum of double- 
Gaussian signal function plus second-degree Chebyshev polynomial background function to 
get the tag mass sideband scaling factor. The invariant mass distributions of tag candidates 
for each tag mode are shown in Fig. [1] and the ST yield and AM(D S ) sideband scaling factor 
are summarized in Table [B We find hst = 26334 ±213 summed over the three tag modes. 

B. Signal-event selection 

A DT event is required to have a ST, a single e + , no additional charged particles, and 
the net charge of the event Q net = 0. We require the momentum of the e + candidate be at 
least 200 MeV. 
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The DT events will contain the sought-after — > r + v T (r + — > e + z/ e ^r) events, but 
also some backgrounds. The most effective variable for separating signal from background 
events is the extra energy (.Eextra) in the event, i.e., the total energy of the rest of the 
event measured in the electromagnetic calorimeter. This quantity is computed using the 
neutral shower energy in the calorimeter, counting all neutral clusters consistent with being 
photons above 30 MeV; these showers must not be associated with any of the ST decay 
tracks or the signal e + . We obtain E extTa in the signal and sideband regions of AM(D S ). 
The sideband-subtracted .Eextra distribution is used to obtain the DT yield. 

The E extra distribution obtained from data is compared to the MC expectation in Fig. [21 
We have used the invariant mass sidebands, defined in Sec. IIII A\ to subtract the combinato- 
rial background. We expect that there will be a large peak between 100 MeV and 200 MeV 
from D* — ► '-/Ds decays (and from D* — ► tt°D s , with a 5.8% branching fraction (3|). Also, 
there will be some events at lower energy when the photon from D* decay escapes detection. 
Based on considerations described in the next paragraph, we define our signal region to be 

Eextra < 400 MeV. 



C. Background estimation 

After the AM(D S ) sideband subtraction, two significant components of background re- 
main. One is from — > K^e + v e decay. If the K° L deposits little or no energy in the 
calorimeter, this decay mode has an E extia distribution very similar to the signal, peaking 
well below 400 MeV. The second source, other D s semielectronic decays, rises smoothly with 
increasing E extra , up to 1 GeV. Estimates of these backgrounds are also shown in Fig. [2j The 
optimal signal region in E extra for DT yield extraction is predicted from an MC simulation 
study. Choosing E extra less than 400 MeV maximizes the signal significance. Note that 
with our chosen requirement of E extra < 400 MeV, we are including — > r + z/ r 7 as signal. 
However, this is expected to be very small, as the kinetic energy of the r + in the Df rest 
frame is only 9.3 MeV and it cannot radiate much. 

The number of nonpeaking background events b np in the E extra signal region is estimated 
from the number of events in the E extra sideband region between 0.6 GeV and 2 GeV, scaled 
by the MC-determined ratio c& (= /b^ 1 ) of the number of background events in the E extTa 
signal region, to the number of events in the E extra sideband region, b^ h \ The number of 
peaking background events b p due to the — > i^°e + z/ e decay is determined by using the 
expected number from MC simulation. The overall expected number of background events 
in the E extra signal region, b, is computed as follows: 

b = b np + b p = c b b {h) (data) + b(K L e + v e ) MC , (7) 

where M h )(data) is the number of data events in the E extTll sideband region and b(Kle + V 6 )mc 
is the number of background events due to —>■ K®e + v P as estimated from our MC 
simulation. We normalize this quantity using our measured [25j B(Df — ► K^e + v e ) = (0.19± 
0.05 ± 0.01)%. We simulate calorimeter response to K° L using a momentum dependent K\ 
interaction probability density function obtained from studying -^(3770) — ► D°D° events 
in which the D° has been reconstructed in hadronic tag modes and the D° decays to the 
K < [tt + tt~ final state. 

The numbers of estimated background events from peaking and nonpeaking sources in 
each tag mode are summarized in Table [Til 



7 



1851208-010 



140 
120 
100 



O 

o 80 



CO 


> 

LU 



60 - 



40 
20 




• Data 

I I MC Total 
I I MC Signal 

MC Background (BG) 



MC D+^K°e + v e 
j MC BG w/o D3 -> K° e + v ( 




• 



• 



1 







0.5 



1.5 



"extra 



(GeV) 



FIG. 2: Distribution of -E ex tra after AM(D S ) sideband subtraction. Filled circles are from data 
and histograms are obtained from MC simulation. The MC signal and peaking background (Df — > 
K®e + i/ e ) components are normalized to our measured branching fractions. The errors shown are 
statistical only. 



IV. RESULTS 

The signal efficiency determined by MC simulation has been weighted by the ST yields 
in each mode as shown in Table IHIl We determine the weighted average signal efficiency 
e = (72.4 ± 0.3)% for the decay chain D+ — > r + u T — > e + v e i> T i> T . 

The DT yields with the 400 MeV extra energy requirement are summarized in Table UVl 
We find noT = 180.6 ± 15.9 summed over all tag modes. Using £>(r + — > e + z/ e z/ T ) = (17.85 ± 
0.05)% 0], we obtain the leptonic decay branching fraction B(Df — > t + v t ) = (5.30±0.47)%, 
where the uncertainty is statistical. 



TABLE II: Estimated backgrounds in the extra energy signal region below 400 MeV in each tag 
mode. Here b p is the peaking background from Df — ► if£e + z^ e decay, b np is the nonpeaking 
background from other D s semileptonic decays, and b is the total number of background events. 
The errors shown are statistical only. 



Tag 


mode 


K P 


K 


b 


DJ 


—* <pn~ 


11.8 ±1.0 


7.7 ±0.5 


19.4 ± 1.1 


Ds 


-» K-K*° 


12.2 ±1.1 


8.7 ±0.5 


20.9 ±1.3 






4.6 ±0.7 


4.5 ±0.3 


9.1 ±0.7 



Total 28.5 ± 1.7 20.9 ± 0.8 49.4 ±1.8 



TABLE III: Summary of the signal efficiency determined by MC simulation. Average efficiency 
e and the tag bias &tag are obtained by using the weighting factor w determined from single-tag 
yields in data. 



Tag mode 


w 


e = cdtAst 


fctag = e'sT/ £ ST 


D~ <M~ 


0.3671 


0.6964 ± 0.0046 


1.0089 ± 0.0058 


Dj -» K-K*° 


0.4154 


0.7337 ± 0.0049 


1.0061 ± 0.0060 


D- -» K-K° s 


0.2175 


0.7536 ± 0.0054 


1.0032 ± 0.0065 


Average 




0.7244 ± 0.0029 


1.0065 ± 0.0036 



V. SYSTEMATIC UNCERTAINTY 

Sources of systematic uncertainties and their effects on the — > t + v t branching fraction 
determination are summarized in Table |V] 

We considered six semileptonic decays, Df — > </>e + z/ e , ^e + z/ e , ?/e + z/ e , K°e + b , e , K*°e + b , e , 
and foe + v e) as the major sources of background in the -E ex tra signal region. The second 
dominates the nonpeaking background, and the fourth (with K^) dominates the peaking 
background. Uncertainty in the signal yield due to nonpeaking background (0.7%) is assessed 
by varying the semileptonic decay branching fractions by the precision with which they are 
known 25|. Imperfect knowledge of B(Df — > K°e + v e ) gives rise to a systematic uncertainty 



in our estimate of the amount of peaking background in the signal region, which has an 



TABLE IV: Summary of double-tag (DT) yields in each tag mode, where n DT is the DT yield in 
the tag mass signal region, n^ T is the yield in the tag mass sideband region, s is the tag mass 
sideband scaling factor, b is the number of estimated background in the extra energy signal region 
after tag mass sideband scaled background subtraction, and «dt is the background subtracted DT 
yield. The errors shown are statistical only. 



Tag mode 


n DT 




s 


b 


n D T 


D~ -> 4>-k- 


79 


1 


0.980 


19.4 ± 1.1 


58.6 ±9.0 


D~ -> K-K*° 


110 


6 


1.000 


20.9 ±1.3 


83.1 ± 10.8 


D~ - K-K° s 


50 


2 


0.999 


9.1 ±0.7 


38.9 ± 7.2 


Total 








49.4 ± 1.8 


180.6 ±15.9 



9 



TABLE V: Summary of sources of systematic uncertainty and their effects on the branching fraction 
measurement. 



Source 


£jneCT on o [ /o ) 


Background (nonpeaking) 


U.7 


DJ — > Kj^e^u e (peaking) 


6.2 


Extra shower 


1.1 


Extra track 


1.1 


<Vnet — u 


i.i 


Non electron 


0.1 


Secondary electron 


0.3 


Number of tag 


0.4 


Tag bias 


0.2 


Tracking 


0.3 


Electron identification 


1.0 


FSR 


1.0 


Total 


4.1 



effect on our branching fraction measurement of 3.2%. 

We study differences in efficiency, data vs MC events, due to the extra energy requirement, 
extra track veto, and Q nct = requirement, by using samples from data and MC events, in 
which both the D~ and satisfy our tag requirements, i.e., "double-tag" events. We then 
apply each of the above-mentioned requirements and compare loss in efficiency of data vs 
MC events. In this way we obtain a correction of 1.6% for the extra energy requirement and 
systematic uncertainties on each of the three requirements of 1.1% (all equal, by chance). 

The non-e + background in the signal e + candidate sample is negligible (0.4%) due to 
the low probability (~ 0.1% per track) that hadrons (tt + or K + ) are misidentified as 
e + [26]. Uncertainty in these backgrounds produces a 0.1% uncertainty in the measure- 



ment of B(D+ — > t + u t ). The secondary e + backgrounds from charge symmetric processes, 
such as 7T° Dalitz decay (n° — > e + e~7) and 7 conversion (7 — > e + e~), are assessed by mea- 
suring the wrong-sign signal electron in events with Q ne t = ±2. The uncertainty in the 
measurement from this source is estimated to be 0.3%. 

Other possible sources of systematic uncertainty include n ST (0.4%), tag bias (0.2%), 
tracking efficiency (0.3%), e ± identification efficiency (1%), and FSR (1%). Combining 
all contributions in quadrature, the total systematic uncertainty in the branching fraction 
measurement is estimated to be 4.1%. 



VI. SUMMARY 

In summary, using the sample of 26 334 tagged decays with the CLEO-c detector 
we obtain the absolute branching fraction of the leptonic decay — > t + v t through r + — > 

e + v e v T 

B{D+ -> t + u t ) = (5.30 ± 0.47 ± 0.22)%, (8) 

where the first uncertainty is statistical and the second is systematic. This result supersedes 
our previous measurement [li| of the same branching fraction, which used a subsample of 
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data used in this work. 

The decay constant fo s can be computed using Eq. t^j with known values 0| Gf 
1.16637(1) x 1(T 5 GeV" 2 , m Da = 1968.49(34) MeV, m T = 1776.84(17) MeV, and r D= 



500(7) x 10~ 15 s. We assume | V cs | = \V ud \ and use the value 0.97418(26) given in Ref. [27 
We obtain 

f Ds = (252.5 ±11.1 ±5.2)MeV. (9) 

Combining with our other determination [15J of fr> a = (263.3 ± 8.2 ± 3.9) MeV with 
— > jj^u^ and — > t + u t (t + — > ti + u t ) decays, we obtain 

f Ds = (259.5 ± 6.6 ± 3.1)MeV. (10) 

This result is derived from absolute branching fractions only and is the most precise deter- 
mination of the D s leptonic decay constant to date. 

Our combined result is larger than the recent LQCD calculation fr> s = (241 ±3) MeV ^ 
by 2.3 standard deviations. The difference between data and LQCD for fn s could be due to 
physics beyond the SM llj], unlikely statistical fluctuations in the experimental measure- 
ments or the LQCD calculation, or systematic uncertainties that are not understood in the 
LQCD calculation or the experimental measurements. 

Combining with our other determination 15j of B(D+ — > t + v t ) = (6.42 ± 0.81 ± 0.18)%, 
via r + — > 7i + u T , we obtain 

B(Dj -> r+z/ T ) = (5.62 ± 0.41 ± 0.16)%. (11) 



Using this with our measurement [15j of B(D+ — > /U + z^) = (0.565 ± 0.045 ± 0.017)%, we 
obtain the branching fraction ratio 

R=^—± —Ll = io.l ±0.9 ±0.3. 12 

This is consistent with 9.76, the value predicted by the SM with lepton universality, as given 
in Eq. ([1]) with known masses 0]. 
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